# cem.output
cases.inweek <- length(cases.in.model$gg10.e)
cases.inweek.matched <- unname(est$tab[1,2])
match.ratio <- unname(est$tab[2,2]/est$tab[1,2])
controltocase.ratio <- unname(est$tab[2,1]/est$tab[2,2])
control.turnout.est <- unname(est$att.model[1,1])
control.turnout.lower <- unname(est$att.model[1,1]-est$att.model[2,1])
control.turnout.upper <- unname(est$att.model[1,1]+est$att.model[2,1])
widower.effect.est <- unname(est$att.model[1,2])
widower.effect.pvalue <- unname(est$att.model[4,2])
widower.effect.upper <- unname(est$att.model[1,2]-est$att.model[2,2])
widower.effect.lower <- unname(est$att.model[1,2]+est$att.model[2,2])
case.turnout.est <- unname(est$att.model[1,1]+est$att.model[1,2])
case.turnout.lower <- unname((est$att.model[1,1]+est$att.model[1,2])-est$att.model[2,2])
case.turnout.upper <- unname((est$att.model[1,1]+est$att.model[1,2])+est$att.model[2,2])

cem.output <- c(cases.inweek.matched, match.ratio, controltocase.ratio, control.turnout.est, control.turnout.lower, control.turnout.upper, widower.effect.est, widower.effect.pvalue, widower.effect.upper, widower.effect.lower, case.turnout.est, case.turnout.lower, case.turnout.upper)

cases.cem.stats <- c(cem.output, cases.in.model.stats)
